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Abstract 

It is now widely accepted that dense filaments of molecular gas are integral to the process of stellar birth and 
potential star-forming cores often appear embedded within these filaments. Although numerical simulations 
have largely succeeded in reproducing filamentary structure in dynamic environments such as in turbulent 
gas and while analytic calculations predict the formation of dense gas filaments via radial collapse, the exact 
process(es) that generate/s such filaments which then form prestellar cores within them, is unclear. In this 
work we therefore study numerically the formation of a dense filament using a relatively simple set-up of a 
uniform-density cylinder in pressure equilibrium with its confining medium. In particular, we examine if its 
propensity to form a dense filament and further, to the formation of prestellar cores within this filament bears 
on the gravitational state of the initial volume of gas. We report a radial collapse leading to the formation 
of a dense filamentary cloud is likely when the initial volume of gas is at least critically stable (characterised 
by the approximate equality between the mass line-density for this volume and its maximum value). Though 
self-gravitating, this volume of gas, however, is not seen to be in free-fall. This post-collapse filament then 
fragments along its length due to the growth of a Jeans-like instability to form prestellar cores like beads on 
a string. We suggest, dense filaments in typical star-forming clouds classified as gravitationally super-critical 
under the assumption of : (i) isothermality when in fact, they are not, and (ii) extended radial profiles as 
against one that is pressure-truncated, thereby causing significant over-estimation of their mass line-density, 
are unlikely to experience gravitational free-fall. The radial density and temperature profile derived for 
this post-collapse filament is consistent with that deduced for typical filamentary clouds mapped in recent 
surveys of nearby star-forming regions. This profile is also in agreement with a Plummer-like density profile. 
For an isothermal filament though, the density profile is much steeper, consistent with the classic density 
profile suggested by Ostriker (1964). On the other hand, increasing the magnitude of the confining pressure 
such that the initial volume of gas is rendered gravitationally sub-critical is unable to collapse radially 
and tends to expand laterally which could possibly explain similar gas filaments found in recent surveys 
of some molecular clouds. Simulations were performed using smoothed particle hydrodynamics (SPH) and 
convergence of results is demonstrated by repeating them at higher resolution. Unlike some of the earlier 
work reported in literature, here we calculate gas temperature by solving the SPH energy equation and 
allow it to cool according to a cooling function. 
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1 INTRODUCTION 

Gas in Giant molecular clouds (GMCs) is distributed 
non-uniformly and appears to aggregate itself into iso¬ 
lated dense clumps or more contiguous and elongated, 
filament-like structures. Detailed observations of poten¬ 
tial star-forming clouds have demonstrated the ubiqui¬ 
tous nature of filamentary clouds. For example, giant 
molecular filaments on scales of a few parsecs have been 
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reported in the inter-arm regions of the Milky-Way (Ra¬ 
gan et al. 2014; Higuchi et al. 2014 and other references 
therein). On the other hand, relatively small (by about 
an order of magnitude compared to the former), dense 
filaments have also been reported within star-forming 
clouds in the local neighbourhood (e.g. Schneider & 
Elmegreen 1979; Nutter et al. 2008; Myers 2009; Andre 
et. al. 2010; Jackson et. al. 2010; Arzoumanian et. al. 
2011; Kainulainen et al. 2011, 2013 and Kirk et. al. 2013 
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are only a few authors among an exhaustive number of 
them). 

Star-forming sites within GMCs are often found lo¬ 
cated within dense filamentary clouds or at the junc¬ 
tions of such clouds and further, these filamentary 
clouds usually show multiplicity, in other words, stri- 
ations roughly orthogonal to the main filament, and 
form hubs (e.g. Palmeirim et al. 2013; Hacar et al. 2013; 
Schneider et al. 2012 & 2010; Myers 2009). In fact, 
inferences drawn from detailed observations of star¬ 
forming clouds have led some authors to suggest that 
turbulence-driven filaments could possibly represent the 
first phase in the episode of stellar-birth, followed by 
gravitational fragmentation of the densest filaments to 
form prestellar cores. Filamentary clouds therefore form 
a crucial part of the star-formation cycle. Consequently, 
a significant observational, theoretical and/or numeri¬ 
cal effort has been directed towards understanding these 
somewhat peculiar clouds. In the last few years we have 
significantly improved our understanding about these 
clouds as a number of them have been studied in dif¬ 
ferent wavebands of the infrared regime of the electro¬ 
magnetic spectrum using sub-millimeter arrays on the 
JCMT and the Herschel (e.g. Nutter & Ward-Thompson 
2007; Andre et al. 2010; Men’schikov et al. 2010; also 
see review by Andre et al. 2014). 

The stability and possible evolution of filamentary 
clouds has also been studied analytically in the past 
and in more recent times. However, these models were 
usually developed under simplifying assumptions. For 
example, Ostriker (1964), developed one of the earliest 
models by approximating a filamentary cloud as an in¬ 
finite self-gravitating cylinder described by a polytropic 
equation of state and derived its density distribution. 
In a later contribution, Bastien (1983), Bastien et al. 
(1991) and Inutsuka & Miyarna (1992) studied the sta¬ 
bility criteria of filamentary clouds under the assump¬ 
tion of isothermality and suggested that such clouds 
were more likely to form via the radial collapse of an 
initial cylindrical distribution of molecular gas. These 
models also demonstrated formation of prestellar cores 
along the dense axial filament via Jeans-fragmentation. 

Formation of dense filaments via interaction between 
turbulent fluid flows has been demonstrated numeri¬ 
cally by a number of authors (e.g. Klessen et al. 2000; 
Bate et al. 2003; Price & Bate 2008, 2009; Federrath 
et al. 2010a; Padoan & Nordlund 2011 and Federrath 
& Klessen 2012). Similarly in other recent contributions 
(e.g. Heitsch et al. (2009); Peters et al. 2012 and Heitsch 
2013), respective authors specifically investigated the 
process that is likely to assemble a dense filament. The 
conclusion of these latter authors supports the idea of 
filament formation via radial collapse of gas followed 
by an accretional phase during which the filament ac¬ 
quires mass even as it continues to self-gravitate. In 
fact, Peters et al. (2012) demonstrated the formation 


of filamentary clouds on the cosmic scale and argued 
that a filament was more likely to collapse radially and 
form stars along its length when confined by pressure 
of relatively small magnitude. In another recent con¬ 
tribution, Smith et al. (2014), have demonstrated the 
formation of dense filaments in turbulent gas, however, 
they have not addressed the other crucial issue about 
the temperature profile of these filaments; observations 
have revealed that gas in the interiors of dense filaments 
is cold at a temperature on the order of 10 K (e.g. Ar¬ 
zoumanian et al. 2010). Heitsch et al. (2013) also ad¬ 
dress the issue of tidal-forces acting on the sides of a 
filamentary cloud having finite length and suggest that 
gas accreted at the ends can possibly give rise to fan-like 
features often found at the ends of infrared dark clouds 
(IRDCs). 

In a semi-analytic calculation, Jog(2013) demon¬ 
strated that tidal-forces also tend to raise the canonical 
Jeans mass, though compressional force-fields are more 
likely to raise the local density and therefore lower the 
Jeans mass. Freundlich et al. (2014) derived a disper¬ 
sion relation for a rotating non-magnetised filamentary 
cloud idealised as a polytropic cylinder with localised 
density perturbations. Under these simplifying assump¬ 
tions, the authors demonstrated that the filament in¬ 
deed developed Jeans-type instability with propensity 
to fragment on the scale of the local Jeans length. These 
conclusions are in fact, consistent with those drawn in 
an earlier work by Pon et al. (2011), or even those by 
Bastien et al. (1991) who arrived at similar conclusions 
from their analytical treatment of the problem. 

On the other hand, Inutsuka & Miyarna (1997), 
showed that a cylindrical distribution of gas is un¬ 
likely to become self-gravitating as long as its mass per 
unit length was less than a certain critical value. Sim¬ 
ilar conclusions were also drawn by Fischera & Mar¬ 
tin (2012) when they performed a stability analysis of 
clouds idealised as isothermal cylinders. In the present 
work we study the dynamical stability of an initially 
non-self gravitating cylindrical distribution of molecu¬ 
lar gas maintained initially in pressure equilibrium. In 
particular, we are interested to see how this distribution 
of gas behaves as it is allowed to cool in the presence 
of self-gravity. We would like to examine if a radial col¬ 
lapse does indeed ensue and the nature of the new state 
of equilibrium attained by the gas. Thus, we would like 
to investigate the physical conditions under which radial 
collapse is likely to be the favoured mode of evolution. 
Although the objective we have set ourselves is similar 
to that of Peters et al. (2012), our work differs on three 
counts : (i) unlike Peters et al. who assumed a polytropic 
equation of state to calculate gas temperature, we cal¬ 
culate gas temperature by solving the energy equation 
coupled with a cooling function, (ii) we do not use a fil¬ 
ament that already is dynamically evolving instead, we 
begin by setting up a cylindrical cloud in pressure equi- 
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librium and then examine if it does indeed collapse to 
form a thin dense filamentary cloud. This relatively sim¬ 
ple set-up represents the intermediate stage of filament- 
formation, in other words, the stage after supersonically 
moving turbulent flows in a molecular cloud (MC) col¬ 
lide to form a filamentary distribution of warm gas (see 
e.g. Klessen et al. 2000; Bate et al. 2003; Price & Bate 
2008,2009; Federrath et al. 2010a; Federrath & Klessen 
2012; Smith et al. 2014), and more importantly, (iii) 
we investigate if gravitational supercriticality is a nec¬ 
essary pre-condition for filament-formation and follow 
the evolution of the dense filament further to study the 
formation of prestellar cores in it. 

We observe that when the magnitude of the confining 
pressure is relatively small such that the initial volume 
of gas is at least critically stable, it does indeed collapse 
radially to assemble a thin dense filament which then 
fragments via a Jeans-like instability on the scale of the 
fastest growing unstable mode to form prestellar cores. 
We also investigate physical properties such as the dis¬ 
tribution of gas density and temperature within the 
post-collapse filament. We demonstrate that the den¬ 
sity profile of the post-collapse isothermal filamentary 
cloud is indeed similar to the classic profile suggested 
by Ostriker (1964); the profile is Plummer-like for a fil¬ 
ament that is allowed to cool dynamically. The article 
is divided as follows : in section 2 we discuss the sta¬ 
bility of a pressure-confined isothermal cylinder. Then 
in section 3 we present the initial conditions used for 
this work and briefly describe the numerical algorithm 
used for the simulations. The results are presented and 
discussed in respectively, sections 4 and 5 before con¬ 
cluding in section 6. 


2 Stability of a pressure confined 
non-magnetic cylinder 

We consider a cylindrical distribution of isothermal 
molecular gas composed of the usual cosmic mixture 
(approximately 90% hydrogen and 10% helium). The 
cylinder has initial radius, R, length, L, and maintained 
at temperature, T gas . In a purely non-magnetic case gas 
within this cylinder is supported by thermal pressure, 
Ptherm , against self-gravity and is confined externally 
by finite pressure exerted by the inter-cloud medium 
(ICM), Picm 0- We assume an initial state of equilibrium 
so that the external pressure, Picm, balances the inter¬ 
nal pressure, Ptherm- The internal and external pressure 
are both thermal by nature since particles representing 
the two media are not imparted any initial momentum. 

The stability of a spherical gas body against self¬ 
gravity is defined in terms of the thermal Jeans mass. 
Now, for a cylindrical distribution of gas its mass line- 


1 By inter-cloud medium we mean here the medium that prevails 
between isolated clumps and filamentary clouds within a GMC. 
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density, M;(r), defined as 

Mi(r = R)= M ( r ^ = j dr 2n rp{r) = ■nR 2 p gas . 

0 ) 

This quantity is analogous to the mass of a spherical 
cloud and as will be seen below, is useful in describing 
the stability of a cylindrical cloud against self-gravity. 
Here, p gas , is the initial density of the cylindrical distri¬ 
bution of gas, assumed to be uniform. A gas cylinder is 
stable against gravitational collapse as long as its mass 
line-density is smaller than the maximum value 

9 „ 2 

{Mi) max = ( 2 ) 

where do is the isothermal sound-speed for the cylinder. 
We define the stability factor, S(r ), as the ratio of the 
mass line-density for a cylinder against its maximum 
value. Thus, 


S(r) = 


M«(r) 


GMi (r) 
2 a 2 0 ' 


( 3 ) 


(Mi(r)) max 

When S(R) exceeds unity, Mi(R) > ( Mi ) max , a radial 
collapse ensues and the cylinder may be described as 
super-critical ; cylinders with S(R) less than unity are 
described as sub-critical (Fiscliera & Martin 2012). We 
describe cylinders with S(R) ~ 1 as critically stable. As 
in our earlier work investigating the stability of star¬ 
less cores (Anathpindika & Di’ Francesco 2013), in the 
following sections we will calculate the stability fac¬ 
tor, S(r ), at different radii within the model cylindrical 
cloud over the course of its evolution to describe its sta¬ 
bility against self-gravity. We now calculate the radius 
of the initial cylindrical distribution of gas. Keeping in 
mind that the pressure exerted by the confining ICM 
balances the thermal pressure within the cylinder, we 
can write 


a 0 Mgas 

ttR 2 L 


so that its radius, 


R = 


/a^Mg 


. 0.5 


V 7 iLP e 


( 4 ) 


Finally, following a simple dimensional analysis 
we define the free-fall time, tff = (GMI -3 ) -1 / 2 ~ 
(Gp) -1 / 2 , which is slightly faster than the free-fall time 
in Jeans analysis; p = p(r), is the density of the post¬ 
collapse filament having radius, r < R 1 and a o, the 
sound-speed of gas within it. A free-fall ensues if tff < 
t sc , where t sc is the sound-crossing time, t sc = T-. Using 
the above inequality we obtain the condition for Jeans 
instability as, 

l > Lj eans ~ (Gp) 1 / 2 ' ^ 

The corresponding mass of the fragment is Mf rag ~ 
MiLjeans, which, using Eqn. (3) above, can be re-cast 
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as 


Atfrag 


2 a 0 P(r = r fi at )S(r = r f i at ) 
(Gp ) 3 / 2 


( 6 ) 


vfiat 5 being the radius of the filament over which gas 
density remains approximately uniform before turning 
over in its wings. 


3 Physical details 
3.1 Initial conditions 


In this demonstrative work we began by placing a non¬ 
self gravitating uniform-density cylindrical distribution 
of molecular gas in pressure equilibrium with its sur¬ 
rounding medium. To this initial distribution of gas we 
assigned a feducial temperature (see Table 1). In this 
work the gas pressure and pressure exerted by the ex¬ 
ternal medium are purely of thermal nature. While am¬ 
bient physical conditions such as the turbulent Mach 
number in a MC could possibly influence filament prop¬ 
erties, we leave this aspect of the problem for future 
investigation. Other free parameters in the set-up were, 
the length of the cylinder, L= 5 pc, contrast, D=10, 
between the initial gas density within the cylinder and 
its surrounding medium, the average initial gas density 
within the cylinder, n gas = 600 cm -3 and the initial 
mass line-density, This choice of gas density 

is also consistent with densities found in a typical star¬ 
forming cloud (e.g. Andre et al. 2013). The temperature 
of the confining medium is calculated by using the con¬ 
dition of pressure-balance at the surface of the cylinder. 
The initial radius of the cylinder was calculated using 
Eqn. (4) above. Though simple, this set of initial condi¬ 
tions is suitable for the extant purposes of studying the 
origin of the density profile and the temperature pro¬ 
file of a typical filament and the formation of prestel- 
lar cores along the length of this filament. Parameters 
used in the realisations performed in this work have 
been listed in Table 1. Sinusoidal density perturbations 
along the length of the initial distribution of gas were 
imposed in five simulations; see Table 1. These pertur¬ 
bations were imposed by revising particle positions, x p , 
to x' p , such that the following expression was satisfied, 


x v = x„ + 


AL 


Jeans 


27T 


sin (kx' p ) 


(Hubber et al. 2006), where the wave-number, k = 
jjjZ- —, the integer, n = 1, 2 for respectively, the funda¬ 
mental mode and the second harmonic; A = 0.1, is the 
amplitude of perturbations. 

The number of gas particles in a simulation are cal¬ 
culated according to - 

Nsas = 8 ' F ' Nneibs ' (/? 7\))’ 


where, Vi n u, is the volume of the initial distribution 
of gas, while the initial average smoothing length, 
Kvg{init), was set equal to 0.25Lj etms /x; X being the 
resolution parameter listed in column 11 of Table 1 and 
defined by Eqn. (10) below; the Jeans length, Lj eans ~ 
0.08 pc, at the density threshold(~ 10 -19 g cm -3 ), 
adopted for representing prestellar cores in this work at 
a temperature of 10 K. Particles representing the initial 
distribution of gas were assembled by extracting a cylin¬ 
der of unit size and length from a pre-settled glass like 
distribution. This unit-cylinder was then stretched to 
the desired dimensions and then placed in a confining 
medium of external-pressure, modelled using the ICM 
particles. The envelope of ICM particles had a thickness, 
30 haygf^nit), along each spatial dimension of the cylin¬ 
der. The ICM particles are special SPH particles that 
exert only hydrodynamic force on the normal gas parti¬ 
cles; gas particles on the other hand, exert both, gravi¬ 
tational and hydrodynamic forces. The cylinder and its 
confining medium were then placed in a box lined with 
boundary particles that prevent SPH particles from es¬ 
caping. Unlike the gas and the ICM particles, the bound¬ 
ary particles are dead and do not contribute to SPH 
forces. The layer of dead boundary particles was there 
to simply hold the gas+ICM particles in the box and 
prevent them from diffusing away. Particles within the 
set-up were initially stationary as no external velocity 
field was introduced. 


3.2 Numerical method : Smoothed Particle 
Hydrodynamics (SPH) 

Simulations were performed using our well tested SPH 
algorithm, SEREN, in its conservative energy mode 
(Hubber et al. 2011). The SPH calculations in this work 
also include contributions due to the artificial conduc¬ 
tivity prescribed by Price (2008), to avoid any spurious 
numerical artefacts possible due to the creation of a 
gap between multi-phase fluids. An SPH particle, in this 
work, had a fixed number of neighbours, N ne n, s = 50. 
Gas cooling was implemented by employing the cooling 
function given by the equation below, 


hjTgas) 

Tc_r 


10 7 exp 


( 


-1.184 x 10 5 \ 

Tgas + 1000 / 


UxlO- 2 ^-^) cm 3 (7) 

(Koyarna & Inutsuka 2002). Constant background heat¬ 
ing due to cosmic rays is provided by the heating func¬ 
tion, Tcfl = 2 x 10 -26 erg s -1 . The temperature of an 
SPH particle is calculated by revising its internal en¬ 
ergy. The equilibrium temperature, T eq , attained by a 
small SPH particle corresponds to equilibrium energy, 
u eq , and the timescale, r, over which energy is either 
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Table 1 Parameters used for each realisation. 
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No. 

[x 10 4 Kcm -3 ] 

[M 0 ] 

[pc] 
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Dynamics 
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1.5 

210 

0.589 

42 

25 

41.37 
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339941 
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1.5 

244.7 

0.636 

48.94 

25 

41.37 

Cooling Func. 
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339941 

1.6 

N 

3 

2.4 

244.7 

0.636 

48.94 

40 

66.19 
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N 

4 
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“Magnitude of the external confining pressure that balances the pressure within the cylinder, Ptherm- 
b Calculated using Eqn. (4) in the text. 

“Initial mass-per unit length of the filament, i.e., the mass line-density. 
d The initial gas temperature. 

“Maximum mass line-density defined by Eqn. (2). 

^Total number of particles including those representing the ICM . 

9 Number of gas particles 
^Resolution factor 

Tf whether perturbations described in §3.1 were imposed in a particular realisation. 


radiated or acquired by a particle is, 

_ U P ~~ ( u p)eq /n\ 

(n 2 A — uTcr) 

Then, after a time-step dt , the energy u of a particle is 
revised according to, 

u = u eq + (u - u eq )exp(-dt/T) (9) 

(Vazquez-Semadeni et al. 2007). However, a more 
accurate calculation of the gas temperature demands 
solving the radiative transfer equation(e.g. Krumholz 
et al. 2007; Price & Bate 2009; Bate 2012), and 
accounting for the molecular chemistry in the cloud 
(Glover et al. 2010). We leave this for a future work. 
The SPH sink particle 

We represent a prestellar core by an SPH sink particle 
in this work, for it is not our intention to follow 
the internal dynamics of a core. An SPH particle is 
replaced with a sink only when the particle satisfies 
the default criteria for sink-formation which include, 
negative divergence of velocity and acceleration of 
the seed particle (see Bate, Bonnell & Price 1995; 
Bate & Burkert 1997 and Federrath et al. (2010b)). 
We set the density threshold for sink-formation at 
Pthresh ~ 10 -19 g cm -3 (~ 10 5 cm -3 ), which despite 
being on the lower side, is consistent with the average 
density of a typical prestellar core. The interested 
reader is referred to Hubber et al. (2011) for other 
technical details related to implementation of this 


prescription in our algorithm SEREN. The radius of a 
sink particle was set at, r S i n k ~ 2.5 h p , where h p is the 
smoothing length of the SPH particle that seeds a sink. 
Sink particles in the simulation therefore have different 
radii. A sink interacts with normal gas particles only 
via gravitational interaction and accretes gas particles 
that come within the pre-defined sink-accretion radius, 

' sink • 


3.3 Resolution 

The average initial smoothing length, h avg ^ n i t ) , for the 
ensemble of calculations discussed in this article is cal¬ 
culated as - 

, _ ( 3 N neibs } !/ 3 

tlavg(iuit) - [7^ ' ' VinitJ , 

where the number of neighbours, Nneibs = 50. The di¬ 
ameter of an SPH particle, d = 4 h avg ( init ) , so that the 
resolution criterion becomes, 

x = (io) 

d 

X > 1, and L jeans, the typical length-scale of fragmen¬ 
tation has been defined by Eqn. (5) above. Evidently, 
higher the value of X, better the resolution. Value of 
the resolution factor, X , for each simulation has been 
listed in column 10 of Table 1. The resolution crite- 
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Figure 1. Cross-sectional rendered density-plots of the mid-plane of the collapsing cylinder in realisation 1. Overlaid on these images 
are velocity vectors representing local gas-motion; shorter vectors indicate relatively lower velocity as compared with those that are 
longer. Time in units of Myr has been marked on the top-right hand corner of each image. 


rion defined by Eqn. (10) above is the SPH equivalent 
of the Truelove criterion (Truelove et al. 1997), defined 
originally for the grid-based Adaptive Mesh Refinement 
algorithm. It has been shown that the gravitational in¬ 
stability is reasonably well-resolved for y = 1, and need 
be ~2, for a good resolution of this instability. Note 
that some authors prefer to define the resolution fac¬ 
tor, A, as the inverse of the fraction on the right-hand 
side of Eqn. (10), in which case the criterion for good 
resolution becomes, X < 1 (e.g. Hubber et al. 2006). 
All our simulations discussed in this work satisfy this 
resolution criterion. 

However, this choice of resolution is reportedly insuf¬ 
ficient to achieve convergence in calculations of crucial 
physical parameters such as the energy and momen¬ 
tum, but sufficient to prevent artificial fragmentation 
of collapsing gas (Federrath et al. 2014). In order to al¬ 
leviate this latter problem these authors suggest a more 
stringent resolution criterion, y > 15. However, our ex¬ 
tant purpose here is limited to studying the stability of 
a filamentary cloud and its fragmentation so that the 
adopted choice of numerical resolution^ ~ 2), should 
be sufficient. In fact, in §4.1 below we demonstrate nu¬ 
merical convergence over the observed fragmentation 
of the post-collapse filament. Finally, the sink radius, 
vsink , defined in §3.2 above is ~ 0.08 pc in simulations 
developed with the least number of particles and ~ 0.05 
pc in Cases 6 and 7, those developed with the high¬ 
est resolution in this work. Although prestellar cores of 
smaller size are known to exist in star-forming clouds 
(see for e.g. Jijina et al. 1999), the resolution adopted 
here is sufficient for a demonstrative calculation of the 
kind presented in this work. The typical width of the 
post-collapse filament, ~0.1 pc, is also consistently re¬ 
solved by a minimum of 80 particles along its length. 

Since this volume of gas is allowed to cool, what fol¬ 
lows, is the fragmentation of cooling gas. The length- 
scale of this fragmentation is comparable to the ther¬ 


mal Jeans length (see for e.g. Vazquez-Semadeni et al. 
2007), which is much larger than the initial width of 
the cylinder so that there is no fragmentation in the ra¬ 
dial direction, as will be demonstrated below. Further¬ 
more, we have set our resolution to the Jeans length, 
Ljeans, evaluated at 10 K, which is lower than the av¬ 
erage gas temperature within interiors of the filament 
as will be evident in the following section. Our simu¬ 
lations therefore have the desired spatial resolution to 
represent core-formation in the post-collapse filament. 

4 Results 

Our interest lies particularly in examining the possibil¬ 
ity of forming a dense filamentary cloud out of a cylin¬ 
drical volume of gas confined by thermal pressure. This 
possibility is explored for different choices of initial con¬ 
ditions; first, when the volume of gas is initially criti¬ 
cally stable, then marginally super-critical and finally, 
sub-critical. Interestingly though, the initial cylindri¬ 
cal distribution of gas evolves radically differently if it 
is at/above critical stability as compared to when it 
is sub-critical. We also vary the spatial resolution and 
compare results from simulations developed by using a 
cooling function with that from an isothermal calcula¬ 
tion. 

4.1 Case 1 (Serial nos. 1, 2, 4, 5, 6 and 7 in 
Table 1; P A"'' — 1.5 X 10 4 K cm -3 ) 

The initial cylindrical distribution of gas, in all these 
test cases, has comparable mass and an identical mag¬ 
nitude of confining pressure, Pi cm . One of these, test 
realisation 1, was developed under the assumption of 
isothermality while the cooling function defined by 
Eqn. (7) above, was used in the remaining test cases. 
The mass line-density, (M ; ) max , defined by Eqn. (2) is 
slightly less than that for the initial distribution of gas, 
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Figure 2. Upper-panel A rendered density plot through the mid¬ 
plane of the post-collapse filament in realisation 1. Black blobs on 
this plot represent the positions of prestellar cores in this filament 
and those that have formed at this epoch are separated by ap¬ 
proximately a Jeans-length for this calculation (~0.2 pc). Lower- 
panel Plot showing a time-sequence of density profile within the 
collapsing filament for this case. The density distribution for the 
post-collapse filament (t = 0.8 Myr), is very well approximated 
by the Ostriker profile. 


(Mi)init, for all the test realisations, except for the first 
where the two are of comparable magnitudes. The initial 
distribution of gas is therefore marginally super-critical 
in all realisations; gas in the first realisation is approxi¬ 
mately in equilibrium in the sense, S(R) ~ 1 (see Table 
1), R being the radius of the cylinder. In all these real¬ 
isations the radius of the initial cylindrical distribution 
of gas begins to shrink due to the onset of an inwardly 
propagating compressional wave. This radial contrac¬ 
tion of the initial cylindrical distribution in outside-in 
fashion is visible in the rendered density images of the 
cross-section through the mid-plane of the collapsing 
cylinder shown in Fig. 1. Overlaid on these plots are 
the velocity vectors that denote the direction of local 
gas-motion. 

This contraction soon assembles a thin dense filament 
along the axis of the cylinder and prestellar cores, rep¬ 
resented by SPH sink particles, begin to form along the 
length of the filament. At this point a word of cau¬ 
tion would be appropriate. By a compressional wave we 
do not imply that the external pressure compresses the 
initial cylindrical distribution of gas in the radial direc¬ 
tion. That hydrostatic cylinders cannot be compressed 
in this manner has been demonstrated by Fiege & Pu- 
dritz (2000). In this work density enhancement along 
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Figure 3. Shown here is the time-variation of the stability factor, 
S(r), defined by Eqn. (3) as a function of the radial coordinate 
within the collapsing filament in realisation 1. Interestingly, gas in 
the collapsing cylinder always remains gravitationally sub-critical; 
S(r) < 1. 



t/[Myrs] 

Figure 4. A plot showing the accretion history of the cores that 
form in the post-collapse filament in the first realisation. 


the cylinder-axis is the result of gas being squeezed dur¬ 
ing the radial collapse of the initial configuration. We 
reserve the attribute of a compressional wave to de¬ 
scribe this process. Shown on the upper-panel of Fig. 2 
is a rendered density plot of the post-collapse filament 
in realisation 1. The black blobs on this image represent 
the location of prestellar cores that begin to form along 
the length of this filament. 

Evidently, the cores that had formed in this filament 
at the time of terminating the calculation are separated 
by approximately a Jeans length(~ 0.2 pc at T gas = 
25 K; avg. density ~ 10 -19 g cm -3 ), defined by Eqn. 
(5). The image on the lower-panel of Fig. 2 is a plot 
of the radial distribution of gas density within the col¬ 
lapsing cylindrical cloud in this realisation. This plot 
and all other similar plots in this paper were generated 
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by taking a transverse section through the mid-plane of 
the filament. The inwardly propagating compressional 
disturbance is evident from the density plots at early 
epochs of the collapse. As the gas is steadily accumu¬ 
lated in the central plane of the cylinder, the density 
of gas there rises and eventually develops a profile that 
matches very well with the one suggested by Ostriker 
(1964), 


p{r) 


Pc 




2-1 p /2 > 


( 11 ) 


with p = 4. Here p c = p(r = rput ); rjiat being the ra¬ 
dius at which the density profile develops a knee. The 
density of this filament falls-off relatively steeply in the 
outer regions which is consistent with that reported by 
Malinen et al. (2012) for typical filamentary clouds in 
the Taurus MC. Furthermore, this figure also shows, 
ffiat ~ 0.1 pc, which is comparable to the radius of typ¬ 
ical filaments found in eight nearby star-forming clouds 
in the Gould-Belt (e.g. Arzoumanian et al. 2011; Andre 
et al. 2014). 

Another interesting conclusion that can be drawn 
from this realisation is that about the stability of a fil¬ 
amentary cloud. It can be seen from Fig. 2 that the 
centrally assembled, post-collapse dense filament at the 
epoch when calculations for this realisation were termi¬ 
nated, has a radius on the order of 0.1 pc, the point 
at which the density-profile develops a knee. An aver¬ 
age density of ~ 10 -19 g cm -3 , the threshold for core¬ 
formation in these calculations, would suggest a mass 
line-density, Mi ~ 50 M 0 pc -1 , for this post-collapse 
filament which is significantly higher than, {Mi(T gas = 
25K)) max ~ 41 M 0 pc -1 . The filament would therefore 
be expected to be in free-fall and should collapse rapidly 
to form a singular line. However, such is not to be the 
case and Fig. 3 shows that the collapsing cylindrical 
cloud, or indeed the post-collapse filament, in fact, re¬ 
mains sub-critical at all radii, r < R. We will discuss 
later in §5.1 the implications of this finding for the clas¬ 
sification of observed filamentary clouds in typical star¬ 
forming regions. 

In Fig. 3 we have shown the radial variation of the sta¬ 
bility factor, S(r), defined by Eqn. (3), and integrated 
over the radius of the collapsing cylindrical cloud. For 
the purpose, using the Sturges criterion for optimal 
bin-size (Sturges 1926), we divided the cylindrical dis¬ 
tribution of gas into a number of concentric cylinders 
having incremental radii in the radially outward direc¬ 
tion. Number of these concentric cylinders is given by, 
N C yi = (1 + ln(IVg OS )); N gas being the number of gas 
particles in a realisation. The innermost cylinder would 
therefore have the smallest radius and radii of successive 
cylinders in the outward direction would be incremented 
by a small increment, dr. Thus, if r, is the radius of the 
i th cylinder, then the radius of the (i + l) th cylinder is 


simply, 'r'i + i = Ti + dr; dr -C R. The initial distribution 
of gas in this case was characterised by S(R) ~ 1. 

Evident from the plot shown in Fig. 3 is the fact 
that the gravitational state of the collapsing cylinder 
does not vary during the process. Within the cylinder 
gas always remains gravitationally sub-critical so that, 
S(r < R) < 1, although it does rise steadily close to the 
centre where the density is the highest in the post¬ 
collapse filament. In other words, shells of gas within 
the collapsing cylinder though self-gravitating, are not 
essentially in free-fall. That this must be the case is also 
evident from the density plots shown in Fig. 1. From 
these plots it can be seen that gas inside the collapsing 
cylinder has a relatively small velocity in comparison 
with that in the outer regions that is pushed in by the 
compressional wave. Also, the observed outside-in type 
of collapse lends greater credence to the conclusion that 
gas within the collapsing cylinder must be gravitation¬ 
ally sub-critical, for had it been super-critical, the col¬ 
lapse would have been of the inside-out type, i.e., the 
innermost regions of the cylinder would have collapsed 
rapidly followed by the outer regions. 

We had presented a similar argument in an earlier 
work (Anathpindika & di’ Francesco 2013), where hy¬ 
drodynamic models developed to explain why some 
prestellar cores failed to form stars despite having 
masses well in excess of their respective Jeans mass, 
were discussed. In that work it was demonstrated that 
cores could indeed contract in the radial direction with¬ 
out becoming singular (i.e., a free-fall collapse was not 
seen in these cores), if the mass of infalling gas at radii 
within the core was smaller than the local Jeans mass, or 
equivalently, we had suggested that the Jeans condition 
had to be satisfied at all radii within a collapsing core 
for it to become singular (protostellar). In the present 
work we have used the mass line-density instead of the 
thermal Jeans mass to discuss the stability of a cylindri¬ 
cal volume of gas. We observe, the collapsing cylinder 
appears to evolve through quasi-equilibrium configura¬ 
tions where thermal pressure provides support against 
self-gravity and a radial free-fall does not ensue. In¬ 
terestingly, we observe that prestellar cores can indeed 
form in a filamentary cloud that is not in radial free-fall. 
These cores acquire mass by accreting gas and shown 
in Fig. 4 is the accretion history of these cores. 

Realisation 2 is a repetition of the previous simulation 
but the gas temperature was now calculated by solving 
the energy equation and further, the energy of a gas 
particle was corrected according to Eqn. (9) to account 
for gas cooling/heating. As with the isothermal cylin¬ 
drical cloud in the previous realisation, the marginally 
super-critical initial distribution of gas in this case also 
begins to collapse in the radial direction. And, as in the 
previous case the inwardly propagating compressional 
wave assembles a dense filament of gas along the axis of 
the initial cylindrical cloud. Shown in Fig. 5 is the radial 
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Figure 5. As in the plot shown in Fig. 1, this is the radial dis¬ 
tribution of gas density within the collapsing filament at different 
epochs in realisation 2. 



Figure 8. Same as the plot in Fig. 3, but now for the cores that 
form in the post-collapse filament in realisation 2. 



Figure 6. As in Fig. 2, shown in this plot is the time-variation of 
the stability factor, S(r), calculated for different radii at various 
epochs of the collapsing filament in realisation 2. 



Figure 7. This plot shows the radial variation of density averaged 
gas temperature at different epochs within the collapsing filament 
in realisation 2. 


distribution of gas density within the collapsing cylin¬ 
drical cloud which, though remarkably similar to that 
shown in Fig. 2 for the isothermal filament, is relatively 
shallow (Plummer-like; p = 2 in Eqn. (11)). The radial 
variation of the stability factor, S(r), within the collaps¬ 
ing cylinder for this realisation has been shown in Fig. 6. 
This plot is qualitatively similar to the one produced for 
simulation 1 and shown earlier in Fig. 3. As with this 
former plot, gas within the collapsing cylinder in this 
case also always remains less than unity, S(r < R) < 1, 
suggesting that gas is not in free-fall. Also, the filament 
in this case is assembled on a timescale comparable to 
that in the isothermal realisation. This is probably be¬ 
cause, despite the gas-cooling, the collapsing gas is still 
gravitationally sub-critical. A direct comparison of the 
mass line-density, Mi(R ), against the threshold for sta¬ 
bility, ( Mi ) max , however, suggests otherwise. This is be¬ 
cause, implicit in a comparison of this kind is the as¬ 
sumption of isothermality and uniform distribution of 
gas density within the post-collapse filament. 

On the contrary, and as is evident from Fig. 7, the 
gas temperature within the post-collapse filament is 
hardly uniform. The cold central region of the filament 
is cocooned by the relatively warm jacket of gas. The 
innermost regions of this filament acquire an average 
temperature of 10 K which is consistent with that re¬ 
ported from observations of filamentary clouds in typ¬ 
ical star-forming regions (see for e.g. Palmeirim et al. 
2013; Arzoumanian et al. 2011; Andre et al. 2010). Fi¬ 
nally, shown in Fig. 8 is the history of core-formation 
in this filament. It is qualitatively similar to that for 
the isothermal filament plotted in Fig. 4. To this end, 
following the formation of the first core, the next set of 
cores form relatively quickly. Also, the filament in the 
second realisation, where gas was allowed to cool, be¬ 
gins to form cores earlier than the isothermal filament 
and in general, cores in the latter are somewhat less 
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the mid-plane of the post-collapse filament (4=1.) in realisation 6 
(only a small volume of the medium confining the filament has 
been shown on this plot). Fine black blobs on top of the image 
represent the positions of cores in the post-collapse filament at 
the time of termination of calculations. Perturbations, in this re¬ 
alisation, were imposed on the length-scale Lj ean3 , defined by 
Eqn. (5). Not all density perturbations have condensed at this 
epoch, but those that are relatively closely spaced are separated 
on a scale on the order of Lj eans ~ 0.08 pc.Lower-panel: Same as 
the picture shown in the upper- panel, but now for the realisation 
4 that was developed with a somewhat lower resolution, though 
sufficient to avoid artificial fragmentation along the length of the 
post-collapse filament. 



Figure 10. Similar to the plot shown in Fig. 8, but with pertur¬ 
bations now imposed on the length-scale, Lj eans /2 (the second 
harmonic), in realisation 7. Again, not all density perturbations 
have condensed at this epoch, but those that are relatively closely 
spaced are separated on a scale on the order of Lj eans / 2 ~ 0.04 
pc. 

massive than those in the former. This is the result of a 
higher mass line-density of the cooling filament relative 
to the isothermal filament in the earlier realisation. 

In realisation 1, for instance, where no external per¬ 
turbations were imposed, we observed that the separa¬ 
tion between cores in the post-collapse filament was on 
the order of the corresponding Jeans length. Next, we 
use the same set of initial conditions as those for reali¬ 
sation 2 in this next set of four realisations namely, 4, 
5, 6 and 7, but now by imposing perturbations on the 
initial distribution of gas as described in §3.1. Perturba¬ 
tions for realisation 4 were imposed on the scale of the 
fragmentation length defined by Eqn. (5), which at 10 
K and average density, ~ 10” 19 g cm -3 , the threshold 
for core-formation, is Lj eans ~ 0.08 pc. We define this 
as the fundamental mode of perturbation. In the next 
realisation, numbered 5 in Table 1, we imposed the sec¬ 
ond harmonic of this perturbation. We repeat this set 
of two calculations with a higher number of particles, 


labelled 6 and 7 respectively, in Table 1. As with the gas 
in realisation 2 that was subject to cooling, in this case 
also, we observe that the initial distribution of gas col¬ 
lapsed radially to form a thin filament aligned with the 
central axis of the cylinder. Shown in Figs. 9 and 10 are 
the rendered images of the post-collapse filaments that 
form in realisations 6 and 7, those that have the best 
resolution in this set of calculations. For comparison 
purposes we have also shown a plot of the post-collapse 
filament from realisation 4 on the lower-panel of Fig. 
9. The physical parameters for this realisation are the 
same as those for realisation 6, but was developed with 
a slightly lower, albeit sufficient resolution to prevent 
artificial fragmentation (Truelove criterion satisfied). A 
comparison of the plots on the two panels of Fig. 9 
suggests little qualitative difference between the out¬ 
comes from the respective realisations. We must, how¬ 
ever, point out that while the satisfaction of the resolu¬ 
tion criterion defined by Eqn. (10) above is sufficient to 
ensure there is no artificial fragmentation, as is indeed 
seen in this work, it is unlikely to be sufficient to obtain 
convergence in calculations of energy and momentum of 
the collapsing gas (Federrath et al. 2014). In fact, be¬ 
low we compare the sink-formation timescales and the 
magnitude of gas-velocity in the radial direction within 
the post-collapse filament. 

Evidently, cores in the post-collapse filament indeed, 
appear like beads on a string. However, we note that not 
all density perturbations had collapsed to form cores 
when calculations were terminated. A few more would 
be expected to form in this filament, but the extant 
purpose of demonstrating a Jeans-type fragmentation 
of this post-collapse filament is served with the set of 
cores that can be seen on either rendered image. Finally, 
shown in Figs. 11 and 12 are the accretion histories 
for cores that form in realisations 6 and 7, respectively. 
There is not much difference between the timescale of 
fragmentation, as reflected by the epoch at which the 
first core appears in either realisation. Although, the 
formation of cores after the first one, in realisation 7 is 
somewhat delayed and form over a timescale of 10 5 yrs, 
as against those in case 6 which form relatively quickly 
after the first core appears. However, it is likely that 
a more significant difference between timescales would 
be visible for a higher harmonic. Also evident from the 
comparative plots shown in Fig. 11 is the absence of 
convergence in the timescale of sink-formation and the 
mass of sink-particles in the two realisations likely due 
to inadequate numerical resolution as was discussed in 
§3.3. Crucially, the number of sink particles remains the 
same irrespective of the choice of resolution which, it 
was argued previously is sufficient to prevent artificial 
fragmentation even for the realisation with the lowest 
realisation in this work. Also, the accretion histories for 
sink-particles in either case are qualitatively similar. 
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Figure 11. Same as the plot in Fig. 8, but now for the cores 
that form in the post-collapse filament in realisation 4 ( upper- 
panel) and 6( lower-panel ). In the former, the realisation with 
relatively lower-resolution, the sink-formation begins a little ear¬ 
lier and sinks are a little more massive than those in the latter 
realisation, one of those that had the highest resolution in this 
work. 



Figure 12. Same as the plot in Fig. 10, but now for the cores 
that form in the post-collapse filament in realisation 7. 



r/[pc] 


Figure 13. The radial variation of the stability factor, S(r), cal¬ 
culated at different epochs for the cylindrical distribution of gas in 
realisation 8. As was seen in Figs. 2 and 5 for realisations grouped 
under Case 1, in this realisation as well the magnitude of S(r ) re¬ 
mains unchanged and gas always remains sub-critical even as the 
initial cylindrical distribution is squashed into a spheroidal glob¬ 
ule. 


4.2 Case 2 (Serial nos. 3 and 8 in Table 1; 

= 2.4x 10 4 K cm -3 ) 
k B ' 

In this second case we raised the gas temperature such 
that the initial cylindrical distribution of gas was ren¬ 
dered sub-critical. We performed one realisation for this 
choice of Pi C m and repeated it with a higher resolution, 
numbered 8 in Table 1. Here we discuss this latter real¬ 
isation. In contrast to the simulations discussed under 
case 1, the radial collapse of the cylinder could not be 
sustained as it became to be squashed by the confining 
pressure into a spheroidal globule that was in approxi¬ 
mate equilibrium with the external medium. The other 
notable feature being, this spheroidal globule was as¬ 
sembled on a relatively longer timescale, t ~ 2.5 Myrs; 
in comparison to that observed in the simulations tested 
under case 1. As with the realisations discussed earlier 
under Case 1, shown in Fig. 13 is the stability factor, 
S(r), for this realisation. This is plot is similar to the 
ones shown earlier in Figs. 3 and 6 in the sense that gas 
within the cylindrical cloud always remains sub-critical 
as S(r < R) < 1. However, it differs from the former 
plots in one respect; the magnitude of S'(r), in the outer 
regions, away from the central axis of the globule, shows 
a significant fall relative to that for the initial configura¬ 
tion. This, as we will demonstrate later, is because gas 
in the outer regions is significantly warmer than that 
observed in realisations grouped under Case 1. 

Shown in Fig. 14 is the density averaged gas temper¬ 
ature within this globule and was made at the time of 
terminating the calculations in this case (t ~ 2.5 Myrs). 
From this plot it is clear that though the interiors of 
this globule are indeed cold, comparable to the cen- 
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Figure 14. Radial distribution of the density averaged gas tem¬ 
perature within the post-collapse globule in realisation 8 at the 
time of termination of calculations (t = 2.5 Myrs). 



Figure 15. A rendered density image showing a projection of 
the mid-plane of the elongated globule in simulation 8. Time in 
the top left-hand corner of the image is marked in Myrs (t = 2.5 
Myrs ). 



Figure 16. Gas within the post-collapse globule shown in Fig. 
14 exhibits a weak velocity field. The inwardly directed gas in the 
outer regions of the globule as it is being squashed has negative 
velocity on this plot. 


tral regions of the dense filament in case 1, gas away 
from the central axis is significantly warmer. The rel¬ 
atively large gas temperature renders it gravitationally 
sub-critical and so, it is unable to collapse radially to as¬ 
semble a thin dense filament along the axis of the initial 
distribution. Instead, it begins to be squashed from ei¬ 
ther side and proceeds to form an elongated spheroidal 
gas-body as can be seen in the rendered density image 
shown in Fig. 15. Overlaid on this image are density con¬ 
tours that readily reveal evidence of sub-fragmentation 
to form smaller cores along the cold central region of 
the globule; velocity vectors on this image indicate the 
direction of local gas-flow within the elongated globule. 

The magnitude of this velocity field has been plot¬ 
ted in Fig. 15. Note that the magnitude of this velocity 
field is consistent with that derived for typical cores 
(e.g. Motte et al. 1998; Jijina et al. 1999). However, the 
local direction of the velocity field is unlikely to reveal 
much about the boundedness of a core which reinforces 
our suggestion made in an earlier work related to mod¬ 
elling starless cores. There it was shown, a core could 
remain starless despite exhibiting inwardly pointed ve¬ 
locity field which usually, is indicative of a gravitational 
collapse (Anathpindika & diFrancesco 2013). Formation 
of sub-structure within this globule can be identified 
with the aid of density contours. This indicates the on¬ 
set of sub-fragmentation. The fragments, at this epoch 
though, have not reached the density threshold for cores 
adopted in this work. 

5 Discussion 

Prestellar cores usually appear embedded within dense 
filamentary clouds. Understanding the morphology of 
these clouds is therefore crucial towards unravelling the 
details of the star-forming process. In the present work 
we explore the possibility of formation of prestellar cores 
along the length of a filamentary cloud. To this end 
we developed hydrodynamic simulations starting with 
a cylindrical distribution of gas and separately studied 
its evolution when it was gravitationally sub-critical and 
super-critical. We identify two possible modes of evolu¬ 
tion : one, when the initial cylindrical distribution of gas 
is either marginally super-critical(5'(i?) > 1), or even 
critically st able (S(i?) ~ 1), a thin dense filament forms 
as a result of radial collapse of the initial distribution 
of gas. The post-collapse filament continues to accrete 
gas during in-fall and prestellar cores form along the 
length of this filament via a Jeans-like instability. In 
literature this mode has been identified as a tendency 
to form a spindle (e.g. Inutsuka & Miyama 1998); and 
two, when the initial distribution of gas is gravitation¬ 
ally sub-critical(S'(/?) <1), it tends to get squashed and 
forms a spheroidal globule that is in approximate pres¬ 
sure equilibrium with its confining medium. The rela¬ 
tively cooler regions of this globule show evidence for 
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further fragmentation into smaller cores. Similar exam¬ 
ples of fragmentation within larger cores to form smaller 
ones has been identified in for instance, the Serpens 
North (Duarte-Cabral et al. 2010), or within a number 
of infrared dark clouds (e.g. Wilcock et al. 2011, 2012). 

Previous numerical work by a number of authors on 
the formation of filamentary clouds, and their stability 
against self-gravity, has culminated in at least three pos¬ 
sible scenarios : (i) Filamentary clouds are often seen in 
magneto-/hydrodynamic simulations of turbulent gas. 
These clouds are believed to be generated due to inter¬ 
action between turbulent gas, followed by enhancement 
of self-gravity once the gas becomes sufficiently dense 
(see for e.g. Klessen & Burkert 2000; Klessen et. al. 
2000; Federrath & Klessen 2013). This model favours 
filament formation via interaction between turbulent 
flows so that there is no global collapse; dense filamen¬ 
tary clouds instead, appear locally (e.g. Federrath et al. 
2010; Peters et al. 2012), (ii) Heitsch et al. (2010, 2013), 
have demonstrated that filamentary clouds could also 
form within a self-gravitating MC, and (iii) fragmen¬ 
tation of gas-sheets is also a possible mode of forma¬ 
tion for elongated clouds. This has been demonstrated 
in both, the simple case of an isothermal sheet (e.g. 
Schmid-Burgk 1967; Myers 2009), as well as in the one 
confined by shocks (e.g. Vazquez-Semadeni et al. 2007; 
Anathpindika 2009; Heitsch 2010). 

Although there is little doubt about the fact that 
dense filaments are a result of fragmentation of larger 
clouds and that prestellar cores usually form within 
such filaments like beads strung on a wire, the pro¬ 
cesses) leading to the formation of a thin filament 
are somewhat unclear and is still a matter of debate. 
Heitsch (2013), for instance has suggested that fila¬ 
ments, during their formation, do exhibit density en¬ 
hancement in their central regions, but this density en¬ 
hancement is not associated with a corresponding re¬ 
duction in the filament-width. Smith et al. (2014), par¬ 
tially agree with this conclusion but have suggested 
that dense filaments are likely to have a static density- 
profile, one that does not vary with time. In other words 
filament-formation is probably a one step process and 
once created, dense filaments remain as they are. In 
simulations grouped under case 1 of this work we have 
demonstrated the formation of a thin dense filamen¬ 
tary cloud via radial collapse of the initial cylindrical 
distribution of molecular gas. Peters et al.( 2012) argue 
that the magnitude of the externally confining pressure 
probably determines if whether a dense filament could 
possibly form. These authors have suggested that a rela¬ 
tively small magnitude of the external pressure is likely 
to support the formation of a dense filament. On the 
contrary, a filamentary distribution of gas is more likely 
to end up as a spheroidal globule when the magnitude 
of confining pressure is relatively large. However, the 
possible cause(s) for such occurrence is(are) not clear 


from their work. Conclusions drawn from the simula¬ 
tions presented in this work are broadly consistent with 
those of Peters et aZ.(2012), but with an added qualifica¬ 
tion about the gravitational state of the distribution of 
gas that precedes a dense filament. This leads us to the 
next question, that about the propensity of a cylindrical 
cloud to collapse radially and assemble a dense filament. 
We observe that when the initial distribution of gas is 
gravitationally sub-critical, apart from a higher magni¬ 
tude of confining pressure, the gas does not collapse in 
the radial direction, instead, it forms a spheroidal glob¬ 
ule. We therefore believe, the gravitational state of the 
initial distribution of gas is likely to hold the key to its 
future evolution. 

We argue, a radial collapse leading to the formation 
of a dense filament is possible when the gas is at least 
critically stable(S(i?) ~ 1), initially. Interestingly, the 
mass line-density for the post-collapse filament in this 
case exceeds its maximum value, ( Mi ) max , for stabil¬ 
ity under the assumption of a uniform temperature of 
10 K, as is indeed seen in the interior of the filament; 
though gas in the wings of the filament is warmer. Un¬ 
der this assumption of isothermality and uniform den¬ 
sity the filament could be described as gravitationally 
super-critical (Fischera & Martin 2012). Although the 
two sets of simulations discussed in this work demon¬ 
strate that an initially sub-critical cylindrical volume of 
gas cannot be induced to collapse in the radial direc¬ 
tion, plots shown in Figs. 2 and 5 suggest that gas in 
a radially collapsing cylindrical cloud does not become 
gravitationally super-critical (S(r < R) < 1). In other- 
words, the mass line-density, Mi, within the cylinder 
never exceeds its maximum value, ( Mi ) max , at any ra¬ 
dius within the cylinder. Simulations grouped under this 
case(Case 1) also demonstrate that a filament need not 
be in radial free-fall in order to form prestellar cores 
along its length. Interestingly, evidence for global in¬ 
fall has recently been reported in the filamentary cloud 
DR21 (Schneider et al. 2010), while Palmeirim et al. 
(2013) and Andre et al. (2014) also report observations 
of some filamentary clouds that exhibit radial collapse. 
Inutsuka & Miyama (1997), by assuming super-critical 
initial conditions and an isothermal gas had demon¬ 
strated such a collapse. However, this latter set of initial 
conditions would always be predisposed to radial col¬ 
lapse leading to the formation of a thin dense filament. 
Despite their propensity to collapse radially, globally 
super-critical filaments are unlikely to be in free-fall. 
The difficulty with the idea of a radially free-falling fil¬ 
ament is that that filament would collapse rapidly to 
form a thin line and it would be a challenge to explain 
the formation of cores in it. It is therefore difficult to 
reconcile the suggestion of free-falling filaments. We will 
revisit this point in the context of observations of fila¬ 
ments in the following subsection. 
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That filamentary clouds are unlikely to experience a 
free-fall collapse can simply be shown by deriving an 
expression for the radial component of gravitational ac¬ 
celeration, g r , within a typical filamentary cloud. For 
the sake of argument, it would be safe to adopt a 
Plummer-density distribution for gas within a filamen¬ 
tary cloud so that the radial distribution of thermal 
pressure within this cloud may be written as, 


dP(r) 

dr 


a o( r ) 


dp(r) 

dr 


+ P(r) 


dao( r ) 

dr 


( 12 ) 


From the Plunnner-density profile we have, 


dp _ —2 p c r 

dr rl (l + (y;) 2 ) 2 ’ 


(13) 


Combining Eqns. (12) and (13), we wind up with an 
expression for the gravitational acceleration as, 


9r(r) 


~ 1 / dP(r) \ 
p(r) V dr ) 


2p c al(r) r _ da%(r) 

P{r)r 2 (1 + (£) 2 ) 2 dr 


(14) 


Then in the limit of r —> r/j 0 t, Eqn. (14) reduces to 


9r 


a o( r = r flat ) > Q 

4 rfi at 


which suggests that filamentary clouds are unlikely to 
have a true super-critical state, i.e., filaments are un¬ 
likely to be in free-fall even if they do collapse radially. 
We have shown that prestellar cores are likely to form 
via a Jeans-type fragmentation of the filamentary cloud 
and have demonstrated this by developing simulations 
with and without initial perturbations to the density 
field. This scenario is consistent with the argument pre¬ 
sented by Freundlich & Jog (2014). We also note that a 
Plummer-like profile (p = 2), fits the radial density dis¬ 
tribution of the post-collapse filament very well when 
gas temperature is calculated by accounting for gas¬ 
cooling. On the other hand, the steeper Ostriker-profile 
(p = 4), fits the radial density profile for the isothermal 
filament. 

Furthermore, the typical radius of the post-collapse 
filament in our simulations at the epoch when calcu¬ 
lations were terminated, is on the order of ~0.1 pc. 
Both these findings are consistent with corresponding 
values derived for filamentary clouds in the Gould-Belt 
(e.g. Arzoumanian et al. 2011; Malinen et al. 2012; 
Palmeirim et al. 2013; Andre et al. 2014). Although a 
few handful number of exceptions have been reported 
with relatively steep density profiles in outer regions 
and have slopes in excess of 2, e.g. the B211/3(p = 2.27) 
filament in the Taurus MC that also has a relatively 
small radius of ~0.04 pc (Malinen et al. 2012; Palmeirim 
et al. 2013). However, it remains to be seen if these post¬ 
collapse dense filaments would continue to self-gravitate 
and acquire even smaller radii; investigation into this 
question is best deferred for a future work. We also note 


that purely hydrodynamic simulations such as the ones 
discussed in this work succeed in reproducing the typ¬ 
ical Plummer-like radial density profile for filamentary 
clouds (also see Smith et al. 2014). This is contrary to 
some of the earlier suggestions about the propensity of 
formation of such filaments in magneto-hydrodynamic 
simulations (e.g. Fiege & Pudritz 2000; Tilley & Pu- 
dritz 2003; Hennebclle 2003). It appears that the slope 
of the radial profiles of filamentary clouds is unlikely to 
be influenced by the presence/absence of the magnetic 
field. 

Also, unlike much of the earlier work we have also 
derived the radial distribution of gas density and tem¬ 
perature for the post-collapse filament and showed that 
is indeed consistent with that derived observationally 
for typical filamentary clouds observed in nearby star- 
forming regions. An interesting point to which atten¬ 
tion must be drawn is about the temporal evolution of 
.the density profile of filamentary clouds. Radial density 
profiles shown in Figs. 2 and 5 demonstrate that the col¬ 
lapse of the initial cylindrical distribution of gas is ac¬ 
companied with a rising central density. This is consis¬ 
tent with the finding of Heitsch (2013). However, we also 
observe that the width of the filament shrinks during the 
process which is inconsistent with the other finding by 
the same author and Smith et al. (2014). As can be seen 
from these plots, the problem could be circumvented 
by adopting a value of the r// a t, the point where the 
density-profile turns-over, smaller than where it is ac¬ 
tually seen to he. Recchi et al. (2014) have argued that 
rotation tends to stabilise a filamentary cloud against 
self-gravity and leads to a density-profile shallower than 
the Ostriker-profile. On the contrary, we have shown 
that even a simple prescription of gas-cooling can in¬ 
deed reproduce a Plummer-like density distribution for 
the post-collapse filamentary cloud. 

On the other hand, upon raising the magnitude of 
external pressure, Pi cm /kB , such that the initial dis¬ 
tribution of gas was rendered sub-critical (because the 
condition of pressure balance at the gas-ICM interface 
demands that gas temperature be raised), we observed, 
such a volume of gas was unable to sustain a collapse 
in the radial direction. Instead, after initially shrink¬ 
ing along the radial direction, gas was squashed lat¬ 
erally and therefore tended to rebound. Confined by 
external pressure, this gas then assembled a pressure- 
supported spheroidal globule on a time-scale larger than 
that for simulations discussed above. Interestingly, Ar¬ 
zoumanian et al. (2013) have recently identified fila¬ 
mentary clouds with a greater internal pressure in the 
IC5146, Aquila and the Polaris molecular clouds. These 
filaments have a relatively large velocity dispersion and 
therefore, are gravitationally unbound. Some of these 
in fact, exhibit signs of lateral expansion with little ev¬ 
idence for core-formation along the filament-axis as is 
usually envisaged in the beads-on-string analogy. On the 
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basis of results obtained from our simulations we sug¬ 
gest, factors controlling the efficiency with which molec¬ 
ular gas cools would determine the propensity to form 
dense filamentary clouds which is simply a restatement 
of the stability argument presented in §2. A larger ther¬ 
mal pressure could have an impact on the molecular 
chemistry of a cloud and therefore, on the efficiency 
with which it could possibly cool. However, in the ab¬ 
sence of a complete treatment of the molecular chem¬ 
istry we cannot address this issue in the present work. 

5.1 Possible implications for observational 
studies of filamentary clouds 

Irrespective of whether the initial cylindrical distribu¬ 
tion of gas ends up as a dense filament, we observe 
that gas in the collapsing cylinder remains gravitation¬ 
ally sub-critical, characterised by the stability factor, 
S(r < R) < 1, at all epochs (see Figs. 3 and 6). This 
suggests, the gravitational state of the gas remains un¬ 
altered over the course of its evolution. Of particular 
interest is the implication of this finding for dense fil¬ 
aments. We have seen that the formation of a dense 
filament is the result of a radial contraction of the ini¬ 
tial cylindrical distribution of gas. This contraction is 
associated with a rise in the central density along the 
axis of the cylinder whence a filamentary cloud is as¬ 
sembled. Although this is our observation in the sim¬ 
ulations discussed here, we note that the radially col¬ 
lapsing gas is not in a free-fall and is therefore char¬ 
acterised by S(r < R) < 1, at all radii within the vol¬ 
ume of gas. This sequence supports the hypothesis that 
filamentary clouds are probably assembled via gravi¬ 
tational contraction during which mass is steadily ac¬ 
creted by the centrally located filamentary cloud. This 
scenario is further reinforced by observational findings 
of accreting filaments found in Taurus and Cygnus X 
(Goldsmith et al. 2008; Nakamura & Li 2008; Schneider 
et al. 2010). 

A comparison of the radial velocity field in the two 
cases discussed in this work would be instructive. For 
Case 1 we take the sixth realisation, one of the two 
in this set with the highest resolution, as a represen¬ 
tative calculation. The radial velocity field in the post¬ 
collapse filamentary cloud has been shown in Fig. 17. 
We note that the magnitude of velocity is relatively 
large in the central region and peters-off in the wings 
of the filament. In fact, the radial distribution of gas 
velocity can be approximated by a power-law of the 
type, V(r) oc r^ 0 4 , which is very close to that for a 
bound object (V(r) oc r~ 0 - 5 ). A similar plot for simula¬ 
tion 8 under Case 2 was shown previously in Fig. 16. 
The difference between these plots is evident; in this lat¬ 
ter case, the gas close to the centre is moving radially 
outward where as that in the outer regions is transonic 
and moving inwards with approximately constant veloc- 



r/[pc) 

Figure 17. Same as the plot shown in Fig. 16 but now for simula¬ 
tions 4 and 6 grouped under Case 1. As usual, a negative velocity 
indicates inwardly moving gas. Although the plots are qualita¬ 
tively similar, the resolution of the simulations developed here is 
insufficient to achieve convergence in the magnitude of velocity of 
the collapsing gas. (See §3.3 for a discussion about the resolution). 


ity (sound speed ~0.55 km/s in this region), as a con¬ 
sequence of being squashed. The findings reported here 
from simulations in Case 1 are consistent with those for 
typical filamentary clouds (see for e.g. Arzoumanian et 
al. 2013). 

Filaments of the type seen in Case 1 are convention¬ 
ally described as gravitationally bound. This brings us 
to the next question - as to whether gravitationally 
bound filaments are also likely to experience a free-fall 
collapse. For the purpose of illustration we consider the 
examples of four filamentary clouds, (i) IC5146 in the 
Cygnus region. This filament with a mass line-density, 
Mi ~ 152 M 0 pc -1 , which is an order of magnitude 
in excess of the maximum mass line-density, ( Mi ) max , 
evaluated at 10 K, and therefore should be expected to 
be in a free-fall collapse in the radial direction. How¬ 
ever, it seems, this is not really the case as the filament 
does not appear to be in radial free-fall (Arzoumanian 
et al. 2011; Fischera & Martin 2012). (ii) The filamen¬ 
tary cloud DR21 is the next example. This is a mas¬ 
sive filament that has two sub-filaments that are grav¬ 
itationally super-critical and show signs of in-fall, but 
no apparent signs of radial free-fall (Schneider et al. 
2010). This filament appears more akin with the out¬ 
come of simulations in Case 2 of this work, (iii) Next, 
we consider the example of the Perseus MC. It has a 
mass line-density, Mi ~ 50 M 0 pc -1 to 100 M 0 pc -1 , 
a good factor of 3-5 larger than the maximum mass 
line-density, ( Mi ) max , for the deduced isothermal tem¬ 
perature of ~12 K (Hatchell et al. 2005). Again, a di¬ 
rect comparison of M; against ( Mi) max for the Perseus 
MC would suggest that the filament is gravitationally 
bound and must collapse radially. However, there is no 
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corroborative evidence to suggest that the filament is in 
free-fall though several cores have been detected along 
its length. In fact, estimates of line-width of the opti¬ 
cally thick C 18 0 emission suggests that that filament is 
probably thermally supported, (iv) Finally, we consider 
the example of the Serpens MC, recent observations of 
which, have revealed evidence of radially infalling gas 
onto the dense filaments in this cloud. Yet, there is no 
evidence to suggest that these filaments are in free-fall, 
however, prestellar cores have been detected along their 
length and appear to have formed on the scale of the lo¬ 
cal Jeans length (Friesen et al. 2014). Results from our 
simulations grouped under Case 1 corroborate these ob¬ 
servations. 

It therefore appears to us that the condition, S(r = 
R) > 1, i.e., (Mi(r=R))^ > 1, is probably necessary to 

induce radial collapse in a cylindrical volume of gas. 
Such a volume of gas could become super-critical by 
accreting mass from its surroundings (Heitsch & Hart¬ 
mann 2014). However, this aspect of the problem is 
beyond the scope of present investigation and there¬ 
fore best left for a future work. Significantly though, 
this collapsing gas is unlikely to be in free-fall and 
S(r < R) < 1 which is consistent with the suggestion of 
McKrea(1957); see also Toci & Galli (2014). In fact, if 
indeed this gas is in radial free-fall, filamentary clouds 
would rapidly end up in a thin line and formation of 
cores in such a collapsing cloud would be difficult to rec¬ 
oncile. On the contrary, a steady in-fall, as gas within 
the collapsing filament is allowed to cool appears to 
be a more promising mode of evolution of filamentary 
clouds where gas is initially at least critically stable. As 
we have seen in the simulations grouped under Case 1, 
such filaments can indeed form cores along their length 
via a Jeans-like fragmentation. This evolutionary sce¬ 
nario also circumvents the uneasy question as to why 
filaments with extremely high column densities, on the 
order of a few times 10 22 cm' 2 , are not found. The likely 
solution perhaps is that filaments are not in radial free- 
fall. Those that are gravitationally bound, shrink in the 
radial direction and eventually acquire a thermally sup¬ 
ported configuration. 

If indeed filaments were to become pressure- 
supported during their evolutionary sequence, it would 
also explain why filaments across star-forming regions 
tend to have comparable widths, on the order of ~0.1 
pc (e.g. Arzoumanian et al. 2011). At 10 K, the temper¬ 
ature typically found in filamentary clouds and an av¬ 
erage density ~10~ 19 g cm' 3 , a representative density 
for prestellar cores, the thermal Jeans length, Lj eans ~ 
0.075 pc, which is comparable to the typical radius of fil¬ 
amentary clouds. We make a similar observation in sim¬ 
ulations grouped under Case 1. For instance, Figs. 2 and 
5, where the radial density profile for the post-collapse 
filament in simulations 1 and 2 have been plotted, 


demonstrate that the filament radius, rput ~ Lj eans , 
at the epoch when calculations were terminated in re¬ 
spective simulations. Although, we note that magnetic 
field could also have a key role to play in the evolution 
of filamentary clouds. The point though, is beyond the 
scope of this article. 

6 Conclusions 

In this work we have demonstrated that an initial cylin¬ 
drical distribution of molecular gas can indeed collapse 
radially to assemble a thin dense filament along its axis 
when the magnitude of confining pressure is relatively 
small such that gas is initially at least critically sta¬ 
ble. In this case a radial collapse ensues and the gas 
is able to cool on a relatively short timescale. Equiva¬ 
lently, it may also be argued that a collapse of this kind 
is likely when the line mass exceeds its critical value re¬ 
quired to maintain stability. However, there is a caveat 
to this argument. This argument is valid only if gas is 
assumed to be isothermal and has uniform density as 
was the case with our initial cylindrical distribution of 
gas. For a typical dense filament this argument cannot 
be applied since its density and temperature distribu¬ 
tion is far from uniform. We argue, a filamentary cloud 
even if gravitationally bound, is unlikely to be in ra¬ 
dial free-fall, instead it is likely to contract in the radial 
direction and consequently, the density along its cen¬ 
tral axis steadily rises as gas is accreted on to it. We 
argue that this evolutionary sequence of a typical fila¬ 
mentary cloud could possibly explain - (i) the dearth of 
filamentary clouds that have extremely high extinction, 
and (ii) why the width of filamentary clouds shows little 
variation across star-forming regions. We suggest that 
filamentary clouds are likely to be thermally-supported 
and their radii(widtlr) are likely to be on the order of 
the local thermal Jeans length (twice the local ther¬ 
mal Jeans length). A Plummer-like density profile ap¬ 
pears to fit very well the radial distribution of gas den¬ 
sity in the post-collapse filament so that magnetic field 
may not be necessary to generate filaments with a rel¬ 
atively shallow slope, as has been suggested in the past 
(e.g. Fiege & Pudritz 2000). We also demonstrate that 
prestellar cores are likely to form via a Jeans-like insta¬ 
bility while it is accreting gas. It also appears, cores are 
likely to form along the length of a dense filament if it is 
at least (gravitationally)critically stable. Processes via 
which filaments acquire mass from their parent cloud 
are therefore likely to hold the key to determining the 
efficiency of star-formation (see also Tafalla & Hacar 
2014) 

On the other hand, when the magnitude of the con¬ 
fining pressure was increased such that the initial distri¬ 
bution of gas was rendered gravitationally sub-critical, 
we observed that a radial collapse was unsustainable. 
Instead, the gas distribution was squashed from either 
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side and demonstrated a propensity to expand laterally. 
The result was an elongated spheroidal globule which 
then showed signs of sub-fragmentation. Such objects 
are more likely to form out of dynamic interactions be¬ 
tween gas flows within molecular clouds. The upshot 
therefore is that potential star-forming dense filaments 
are likely to be found in molecular clouds where am¬ 
bient conditions contrive to render the volume of gas 
immediately preceding the dense filament critically sta¬ 
ble, at least (i.e. mass line-density approximately equal 
to its maximum stable value). The possible effect of the 
ambient conditions on the chemistry within a molecu¬ 
lar cloud must directly impact the gas thermodynamics 
and therefore, the dynamical stability of that cloud. We 
will examine this issue in a future contribution. 
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